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GALACTIC AND EXTRAGALACTIC STUDIES, XXI._ DISTRIBU- 
TION OF 95,000 GALAXIES IN THE NORTH GALACTIC CAP 


By HARLOW SHAPLEY 
HARVARD COLLEGE OBSERVATORY 
Communicated March 1, 1951 
Among the regions in which surveys of galaxies have been made in the 
course of a census of the inner Metagalaxy is the important North Galactic 
Polar Cap, covering the 2760 square degrees from 8 = +60° to 8 = +90°. 
The distribution in this area, heretofore unpublished, is now summarized 
for galaxies brighter than photographic magnitude 17.5. The deeper sur- 


veys, now in progress with Schmidt cameras, should soon provide sampling 
counts to at least the nineteenth magnitude. 


1. The galaxy counts summarized in table | were made on plates of 
two series, all of three-hour exposures. The 57 A plates, made with the 
Bruce telescope in South Africa, show stars to magnitude 18.2 on the aver- 
age; the star limit for the 69 MC plates, made at the Oak Ridge Station 
with the Metcalf doublet, averages 17.9. In declination the Polar Cap 
extends from —2° to +58°. North of declination +22° the MC instru- 
ment was used and the area south of +22° (about 40 per cent of the Polar 
Cap), except for MC 33360 at + 19°4, was covered with the southern 
telescope. 

Earlier investigations with these instruments show that the galaxy 
“completeness” limit m, is 0.4 mag. brighter than the star limit m, for 
the A plates,' and 0.57 mag. brighter for MC plates.?. The difference m, — 
m, varies from plate to plate because of slightly varying conditions in seeing 
and focus. The survey therefore lacks uniformity because of this spread 
and also because of the residual uncertainties in the magnitude sequences 
which are based as usual on the Seares and van Rhijn star-count tables. 

2. In order to obtain a uniform picture of the density distribution of 
galaxies to magnitude 17.5, it has been necessary to reduce the galaxy 
counts on each plate by the application of the factor f = 109° 
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PLATE 
NUMBER 
A6718 
A6719 
A6720 
A13810 
Al4551 
Al4609 
Al5251 
A15363 
A15370 
A15384 
A16086 
A16090 
A16407 
A16757 
A1I7345 
A17468 
Al7721 
A17761 
AI7774 
A18201 
A18296 
A18391 
A18506 
A19063 
Al19349 
Alg422 
A19522 
A200905 
A20107 
A20115 
A20117 
A20670 
A20705 
A20862 
A20880 
A20890 
A20918 
A20022 
A21710 
A21926 
A22393 
A22464 
A22560 
A22575 
A23089 
A238153 
A24823 
A26301 
A26304 
A26704 
A26857 
MC27104 
MC27134 
MC27509 
MC27639 
MC 27641 
MC27688 


TABLE 1. 
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GALAXIES IN NortH GaLactic Cap 
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PLATE 
NUMBER 
MC27704 
MC28714 
MC28754 
MC28788 
MC 28868 
MC28914 
MC28938 
MC29265 
MC29283 
MC20316 
MC20377 
MC29382 
MC20447 
MC20404 
MC20408 
MC29526 
MC20530 
MC29532 
MC20606 
MC29612 
MC29945 
MC29949 
MC29962 
MC20078 
MC29979 
MC29981 
MC20986 
MC29987 
MC29988 
MC299080 
MC30038 
MC30041 
MC30045 
MC30129 
MC30734 
MC30740 
MC30797 
MC30837 
MC308F 1 
MC30879 
MC31292 
MC31358 
MC31350 
MC31365 
MC31368 
MC31431 
MC31432 
MC31440 
MC31442 
MC31444 
MC32130 
MC32148 
MC32157 
MC32226 
MC32281 
MC32284 
MC32867 
MC32872 
MC32908 
MC32914 
MC33347 
MC33360 
MC33924 


GALACTIC 


LONG. 
98° 
106 
115 
129 
120 
63 
100 


LAT 
+74 
+57 
+67 
+68 
+58 
+69 
+ 58 
+64 
+80 
+68 


— 


— w& 


| 
No ms Nr No ms Nr 
245 285 17.8 36.4 206 18.1. 37:0 
258 +80 366 17.8 468 319 18.1 34.0 
262 75: 204. 17.9 28:7 339 17.8 64.3 
255 +70 521 179 S579 216 18.1 23.0 
232 +77 423 17.9 47.0 185 18.3 15.0 : 
257 +04 173 17.9 19.2 194 17.8 31.2 
251 +64 0617.1 «31.9 279 18.0 34.1 
280 +80 88 17.1 296 142 420 18.4 29.4 
; 278 +74 125 176 21.0 137 273 18.5 16.7 
296 +73 158 17.7 23.3 198 568 18.1 60.6 
249 +69 278 181 23.4 70 161 17.8 25.9 
: 263 +60 411 17.5 795 185 +69 279 18.1 29.8 
344 +62 241 17.9 268 58 +64 179 17.9 25.1 
341 +69 347 180 33.6 10 +70 357 17.7 66.2 
324 +59 231 17.8 296 73 +57 283 17.8 45.6 } 
267 +65 726 180 703 303 +85 283 18.0 34.6 
194 +66 329 17.5 63.6 354 +69 437 17.9 61.2 
335 +74 500 179 556 341 +80 454 18 L 48.4 
298 +70 1168 18 2 857 52 +50 163 17.8 262 
328 +60 362 18.3 23.1 32 +65 304 17.3 98.0 
313 +60 604 184 33.6 140 +61 288 18.3 23.3 
- : 311 +68 503 18.5 24 138 +71 251 18.0 30.7 
255 +50 441 18.2 32 47 +65 233 18.3 18.6 
180 +62 312 18.0 30 146 +66 306 18.0 37 2 
336 +64 538 186 22 41 +61 142 18.1 15 
212 +60 424 17.9 i7 of 27 +74 464 17.8 74 
218 +78 862 17.8 110 124 +76 390 17.8 62 
t 320 +71 670 184 37. 31 +59 268 17.7 49 
286 +72 876 180 84 158 +62 186 17.8 30 
204 +65 668 18.3 42 22 +68 #198 17.5 48 
321 +78 540 18.0 52 74 +74 447 180 51 
304 +63 618 182 45 83 +59 210 17.7 39 
210 +67 628 18 1 53 170 +62 170 17.8 27 
351 +62 387 176 64 55 +70 320 18.1 35 
220 +71 673 17.6 112 73 79 288 18.0 35 
281 +66 657 184 36 131 +58 612 18.1 65 
226 +60 548 18.3 35 97 +64 366 17.8 59 
198 +70 918 184 51.0 72 +63 169 17 6 35 
219 +64 582 18.1 49.2 85 +64 394 17.8 63 
269 +70 1635 18 6 69.0 102 +79 145 18.0 17 
276 +65 1092 18.8 34.9 125 +63 192 17.6 40 
268 +60 796 18.6 33.6 230 124 17.7 23 
241 +62 658 18.3 42.0 85 +69 212 17.8 34 
: 320 +64 1766 18.6 74.5 158 +86 344 18.0 42 
206 +81 1435 18.5 69.7 60 +83 124 17.7 23 
203 +63 882 18.3 56.2 166 +81 299 17.5 73 
227 +67 834 18.7 30.6 142 +79 360 17.7 66 
206 +74 944 18.4 52.4 10 +80 176 17.9 24 j 
206 +63 278 18.4 15.4 35 +79 413 17.7 76 
; 280 +64 418 18.4 23.2 50 +75 214 17.5 52 i 
306 +74 240 18.2 17.6 154 +74 411 17.6 87 
327 +68 861 18 0 83.0 174 +76 196 17.5 47.9 
285 +62 488 18.5 23.7 159 «+66 224 17.6 47.6 : 
| 332 +57 548 18.5 26.7 358 +65 471 17.9 65.9 
' 230 +74 1712 18.8 54.8 i) +68 173 17.6 36.7 
340 +58 622 184 34.6 191 +75 277 17.4 77.6 
238 +70 «130 18.0 416 172 +68 564 17.8 909 
105 +65 336 18 4 23.6 180 +63 497 17.7 92.2 
16 +87 838 18.5 51.2 350 +75 462 17.8 74.4 
\ 111 +84 277 18.1 29.6 10 +63 363 17.9 508 
4 1145 +72 371 18.3 30.1 4 +75 542 18.0 66.2 i 
108 +62 614 18.5 37.6 2u4 +81 166 18 0 203 
95 +59 262 17.9 36.7 20 +62 250 17.4 70.0 
192 
{ 
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for the MC plates, and f = 10°-°'7-’ ~ "© for the A plates, where m,, in the 
fifth column of table 1, is the observed stellar magnitude limit. 

It follows that Ny, = fNs, where Vy is the total number of galaxies found 
on each plate in a central area 3° on a side, and Ng, is the number we would 
expect to find if the respective plates gave a complete census of galaxies to 
magnitude 17.5, and no fainter. Applying these reductions to Ng in the 
fourth column of table |, and dividing by 9, we get the number of galaxies 
per square degree ( N, in the last column), reduced to magnitude 17.5. 


FIGURE 1 
For example, when an A plate shows stars to magnitude 18.3, and 362 
galaxies are counted in the central 9 square degrees, we find 
No, = 362-10°-607.9 18.3) 908, 


and N, = 23.1 is the computed average number of galaxies per square de- 
gree for that plate. 


4 
90° = 
Om 
270° 
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In making this reduction we have assumed, as heretofore, that we have, 
on the average, a uniform radial frequency of galaxies. There is, as figure 
| shows, marked unevenness in the frequency from plate to plate. Fre- 
quently there is irregular distribution on the individual plates. The 
method is therefore dependable only for the mean of many plates. It is the 
best we can do in the absence of magnitudes for the individual galaxies. 

On 11 MC plates the magnitudes of galaxies 17.5 and brighter have been 
directly determined in earlier studies.* The difference AN between the 
measured galaxies per square degree for these MC plates and the com- 
puted numbers of table 1 for the same plates, averages 3 per cent, or less 
than one galaxy per square degree, the computed numbers being slightly 
larger. 


TABLE 2 
REGIONAL DISTRIBUTION 
NUMBER 
OF 
REGION rorTaL Nor PLATES 
AN’ to 11,468 
AMO to 11,214 
to A270 15,324 
A270 to A360 13,140 
AO to 22,682 
to 28,464 
490 to A270 26,538 
4270 to 24,608 
8 +57 to 8 +70 27,855 f o 
8 +70 to B +80 18,103 4 
B +80 to B +90 5,188 é om 


All plates 51,146 


3. For figure | we have indicated the essential completeness of the sur- 
vey over the cap by taking the area for each plate equivalent to 25 square 
degrees. The plates actually cover nearly 35 square degrees, but for the A 
plates the quality for nebular discovery and measurement falls off outside 
the central 9 square degrees, and for the MC plates outside the central 25 
square degrees. The MC plates require essentially no correction in their 
central 25 square degrees.’ If we had used these larger MC areas, the 
values of .Y, would not have been changed appreciably except for the 
plates: 

25 sq. deg. 9 sq. deg. 
MC27134 where N,; would be 33.6 instead of 51.2 
31368 23.0 
31440 24.6 
32148 47.9 
32872 34. 92.2 


ere 
my 
14.73 
14.80 
14.66 
14.77 
14.76 
14.71 
14.72 
14.75 
14.80 
14.60 
14.80 
14.74 
| 
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In the mean, for all 69 MC plates, we have 
Nr — _ _ gy 

From this small percentage difference we conclude again that distance 
correction on MC plates can be ignored in a galaxy census, and that the 
total number of galaxies entering the survey of the Polar Cap can appro- 
priately be increased from the 55,686 of column 4, table 1, to 94,880. In 
addition, the outer unused parts of the A plates show approximately $5,000 
galaxies, mostly brighter than magnitude 17.0. 

4. Figure 1 shows the irregularities in surface distribution that are now 
generally recognized as characteristic of the inner Metagalaxy. When ex 
tensive areas are considered the unevenness is largely smoothed out. Table 
2 illustrates the comparability of the four quadrants of galactic longitude, 
and the still closer similarity when the Cap is divided into halves. When 
the means are taken for intervals of galactic latitude, the greatest richness 
is found between latitudes +70° and +S0°, as figure | also indicates. 

Scattered over the Polar Cap are several clusters of galaxies, mostly with 
their individual members fainter than magnitude 17.5. A few of these 
were known before the present survey. No attention is paid to them at 
this time because the Schmidt cameras can soon outline and analyze them 
far better after adding two magnitudes to the depth of the survey. 

Finally it should be noted that this North Galactic Cap is much richer 
to magnitude 17.5 than we have found (Harvard, Mount Wilson, Lick 
surveys) for the high-latitude sky at large. In the formula for uniform 
radial space density, log NV = 0.6(m — m,), the space-density parameter, 
m,, for the whole sky falls between 15.1 and 15.2. (This is the magnitude 
to which the average survey must reach in order to find one galaxy per 
square degree.) The last column of table 2 shows that in the North Galac- 
tic Cap m, is 14.74, that is, four-tenths of a magnitude brighter, corre- 
sponding to 1.75 times the average density. In the South Galactic Polar 
Cap we found some years ago' m, = 15.16 * 0.02 (m.e.) from a study of 
some 13,000 galaxies. This important inequality in registered population 
between north and south apparently cannot be attributed to differences in 
space absorption near the galactic plane (with the south having much 
heavier absorption). When galaxies to magnitude 19.0 are examined, the 
inequality between the north and the south has largely disappeared, and 
with it the hypothesis that local space absorption is responsible. The 
inequality was pointed out earlier‘ as one of the indications of the existence 
of significant radial and transverse density gradients in the distribution of 
faint galaxies. 

Of the several assistants who have worked on the survey of the North 
Galactic Cap, Miss Ann Hearn has contributed most in bringing the investi- 
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gation toa conclusion. The MC plates were all made by Mr. H. A. Sawyer, 
and the Bruce plates under the direction of Dr. J. S. Paraskevopoulos. 

Summary. Working with long-exposure plates made with the Bruce 
and Metcalf doublets at the South African and Oak Ridge Stations of the 
Harvard Observatory, we have investigated, in the course of the survey of 
the inner Metagalaxy, the distribution of approximately 95,000 galaxies 
brighter than magnitude 17.5 in the North Galactic Polar Cap. 

About 95 per cent of the area with galactic latitudes greater than +60° 
is covered by 126 plates, but in only 75 per cent of the area have the galaxies 
been tabulated, since the counts could not be used dependably for the 
outer portions of the plates. 

The average population of galaxies over the whole Cap, 45.1 per square 
degree, to magnitude 17.5 (which corresponds to a limiting distance in the 
polar caps of about 80 million light years for the galaxy of average luminos- 
ity) appears to be about 1.75 times that of the South Galactic Cap. The 
distribution diagram and tables | and 2 present fully the survey and the 
irregularities. 

The important Virgo and Coma clusters of galaxies lie within the area 
surveyed, but the latter is too faint and the former too sparse to make 
significant contributions to the census. 


' These PROCEEDINGS, 26, 169 (1940), Harvard Reprint 194. 

2 Harvard Annals, 106, 4 (1938); these PROCEEDINGS, 26, 560 (1940), Harvard Re- 
print 208. A mean value has been used in this paper. 

‘ Seyfert, C. K., Harvard Annals, 105, No. 10 (1937); these PROCEEDINGS, 26, 554- 
561 (1940), Harvard Reprint 208. 

* American Scientist, 32, 67 (1944), Harvard Reprint 255. 


CURVE OF GROWTH POLARIZATION 
By JAMES W. WARWICK 
Wuitrn OBSERVATORY, WELLESLEY COLLEGE * 
Communicated by D. Menzel, February 28, 1951 


Introduction..-When gas is submerged in a magnetic field, the atomic ab- 
sorption lines split into components, the amount of the splitting being pro- 
portional to the strength of the magnetic field. The conditions of polari- 
zation of the absorbed light vary according to the angle between the mag- 
netic field and the direction of propagation of the light. For each com- 
ponent, one can imagine the incident light split into two orthogonal states 
of polarization, one of which cannot be absorbed by the component. 

Magnetic fields of the strength observed in sunspots and stars are fre- 
quently strong enough to separate completely the components of normal 
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triplets lying on the 45° or flat parts of the curve of growth. Each com- 
ponent then acts independently of the others on its own two states of 
polarization, passing one completely, and absorbing only the other. 

In case of complete separation, one can form the curve of growth sepa- 
rately for each component, if he knows the profile of the component as a 
function of frequency and effective numbers of atoms. By means of a 
suitable polarizing analyzer one can eliminate, for any given direction of 
observation, the state of polarization which is passed by a given component. 
One can then determine the profile of the component in units of the con- 
tinuum of polarized light in the state of polarization which can be absorbed 
by the component. 

We shall assume that the equivalent 
width of the component, measured in 
this way in a polarized continuum, de- 
pends on the absorption coefficient in 
the component in a simple, idealized 
manner. The absorption coefficient, in 
turn, depends on the angle made by 
the line of sight with the magnetic field. 
For all angles such that the optical 
depth at component center, X, is less 
than a certain value, X», we assume 
that the equivalent width, W, varies 
proportionally toX. When X is greater 
than Xo, we assume that W is con- 
stant, and equals Wo (Fig. 1). Assumed curve of growth. 

Babcock! used this idealized curve of 
growth in his analysis of the magnetic intensification of stellar absorption 
lines. The justification for this procedure is that the ordinary curve of 
growth’ for unpolarized lines varies linearly with optical depth for small 
enough numbers of absorbing atoms, and comes to a point of inflection for 
moderate numbers of absorbing atoms. The slope at the point of inflection 


Xx 


FIGURE | 


is small. 

It is not necessarily true that the equivalent width of a component, con- 
sidered as a function of its directionally dependent absorption coefficient, 
will vary in the same way as an unpolarized line, considered as a function 
of optical depth. However, in either case, the errors made in assuming a 
simplified form of the curve of growth are probably greater than the error 
of considering polarized components and unpolarized lines to have curves 
of growth of the same shape. 

We shall assume conditions of saturation such that the ¢ components of a 
triplet lie, for all directions, on the flat part of the curve of growth. This 
means (Fig. 1) that each ¢ component has an optical depth greater than Xo, 


> 
| | 

i| 

| 
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for any direction of observation. On the other hand, the single  compo- 
nent of the triplet does not absorb at all in the longitudinal effect. The x 
component reaches the point Xo for some direction (determined by the 
saturation of the o components) intermediate between the longitudinal and 
transverse directions. From that direction up to the transverse direction 
its optical depth increases, but its equivalent width remains constant at W). 

In the case of the triplet, there are two 7 components seen transversally 
and one m component. All three components have the same equivalent 
width, I, in our special conditions of saturation. Each o component ab- 
sorbs light from the polarized continuum characterized by an electric vector 
vibrating in the direction perpendicular to the magnetic field. The 7 com- 
ponent, on the other hand, absorbs from the orthogonal continuum, char- 
acterized by an electric vector vibrating in the direction parallel to the 
magnetic field. . 

If we block out the continuum parallel to the magnetic field by means of 
a Nicol prism, we may expect to observe an equivalent width 2W,. If we 
block the continuum perpendicular to the field, we may observe equivalent 
width I. ‘The equivalent width of the line, as measured through a Nicol 
prism, depends on the orientation of the analyzer. This effect is called 
“curve of growth polarization.” 

Computation of Stellar Equivalent Widths.-The optical depth in either 
o component is given by X, (1 + cos* y), where y is the angle between the 
line of sight and the magnetic field. X, is the optical depth in either ¢ 
component, when y = 2/2. The assumed conditions of saturation are 
such that X,; > Xo. 

If the conditions of saturation had been chosen so that the # and both o 
components had optical depths less than Xo, for the transverse case, then 
no curve of growth polarization would have been present. The optical 
depth in the r component varies as sin® y. For the present case of satu- 
rated components we therefore write the optical depth in the r component 
as 2X, sin® y. 

The optical depth in the x component is X» when sin? y = Xo/2X). The 
choice of a value of Vo, X, determines the saturation conditions. For rea- 
sons of convenience, we take XYy/X; = '/». Then, for sin y < !/s, the x 
component lies on the part of the curve of growth for which W increases 
proportionally to X. 

On the proportional part of the curve of growth, WoW) = X/Xo. If 
we define the unit of equivalent width to be Wo, and the unit of optical 
depth, Xo, we have W = X. The unit of equivalent width is measured in 
the continuum of polarized light which can be absorbed by the component. 
Therefore, the equivalent widths in figure | must be divided by 2 to refer to 
the total, unpolarized continuum. 

The angle 4 is defined to be the angle made between the projection of the 


| 
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magnetic field on the plane of the sky, and the direction of the magnetic 
dipole axis, which we shall assume lies in the plane of the sky. Following 
Babcock, we write the stellar magnetic field, as a function of position on the 
surface of the star, in the form 
IT, = (3/2) H, sin @ cos @ cos ¢ 
H, = (8/2) H, sin @ cos @ sin ¢ (1) 
H, = H,(cos? 6 — '/2 sin? 6). 
IT,,y,2 are the Cartesian components of the magnetic field; //, is along 
the dipole axis and //, along the line of sight. @ and ¢ are spherical coordi- 
nates defined in the usual way with respect to the Cartesian coordinates. 
We introduce the parameters 
10 + 6 cos 26 
a= (2) 
(3 sin 
and 
(1 +- 3 cos 26)? 
(3 sin 20)? 


in terms of which one can verify that 


cos® ¥ 


(5) 


Also, 
(6) 


We have 


y cos? 6 = (7) 


The stellar radiation field is assumed to be of the form 
I =c(A + By), (8) 
where we set the radiation from the center of the disk equal to ¢ by taking 
A+B= 1. (9) 


We normalize J so that the total flux from the visible hemisphere equals 
unity. Thus, we determine ¢ so that 


+ Buju du = 1, 


3 l 
(3 — B) 


cos’ ¢ 
| - (4) 
a 
and 
sin? 
tan? = = 
B 
or 
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One may write the equivalent widths of the coraponents, in terms of the 
angles y and 4, in the forms presented in table 1. The entries of table | 
are the equivalent widths of a single component as measured in the un- 
polarized continuum in the state of polarization parallel or perpendicular 
to the magnetic dipole axis. 

TABLE 1 
LocaL EQUIVALENT WIDTHS OF COMPONENTS 
PAKALLEL PERPENDICULAR Y 


1 
component cos? 6 = sin? 6 y> r/6 


component 2 sin? y cos? 6 2 sin? y sin? 6 
1 cos? y cos? 6 + sin? 6 1 cos? y sin? 6 + cos? 6 


all y 
2 1 + cos? y 2 1 + cos? y 


component 

The longitudinal effect of the o components represents the equivalent 
width the line would have had had there been no Zeeman splitting at all, 
since the two o components absorb in orthogonal states of polarization in 
the longitudinal effect. Because of the flux normalization of the radiation 
field, the integrated equivalent width over the whole stellar disk would 
have been unity, if there had been no Zeeman splitting. The appropri- 
ately weighted sums of the equivalent widths in table | will therefore repre- 
sent the intensification factor introduced by Babcock. 

Figure 2 represents the projected stellar disk. 
Within the hatched areas, y < 7/6. For @ lying 
in the interval 6; to 6, the integration of the 
expressions of table 1 is in two parts, corres- 
ponding to y 2 9/6. 

There are two practical ways to find the nu- 
merical values of the integrals arising from table 
1: (1) cubature, or (2) analytic integration 
over y, followed by quadrature in @. I have 
followed the latter course. The analytic inte- 
grals are all elementary, but somewhat tedious 
and lengthy, and will not be reproduced here. 
The functions resulting from these integrations 
were evaluated for twelve equally spaced in- 
tervals in each of the ranges 0° to 6), 6 to # and @ to 90°. The @ cor- 
responding to the maximum lateral extent of the cross-hatched area rep- 
resents a latitude of are sin 43/3 = 35°3. This is the latitude at which 
the meridional magnetic field is pointed at the observer. The projected 
height of the zone between 6, and 6 is 0.416 radii, and its projected width, 
1.000 radii. 

The equivalent widths were obtained from Weddle's six-interval rule 


FIGURE 2 


The projected stellar disk. 
Regions for which y < 1/6 
are hatched. 


| 
| 
4 ; 
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of quadratures. The formulae expressing the equivalent widths measured 
in the perpendicular and parallel states were found to be 


W, = 


1 
(3 —B) (2.19 — 0.65B) (11) 


1 
W, = B (2.07 — 0.71B). (12) 


(3 
Typical values of B are used to form table 2, whose last column contains 
essentially the polarization. 

Summary and Conclustons.—It appears that 4 or 5 per cent plane polari- 
zation is quite likely to be present in the absorption lines of the magnetic 
stars for some directions of observation. This prediction is based on: 
(1) an idealized curve of growth theory; (2) well-separated Zeeman com- 
ponents; (3) the stellar magnetic field oriented like that of a dipole; (4) 
the star a sphere; (5) darkening to the limb following the cosine law; (6) 
the root-mean-square atomic velocity being constant over the surface of the 
star; (7) the magnetic dipole axis being in the plane of the sky; and (8) 
the line being a Zeeman triplet. 

TABLE 2 
POLARIZATION OF A TRIPLET LN A MAGNETIC STAR 
B Wi + Wi - Wi, 
0.00 1.42 0.04 


0.50 1.48 0.06 
1.00 1.45 0.09 


Whether such a small amount of polarization is actually observable is 
difficult to say. Struve, in attempting to observe resonance polarization 
in the ionized calcium emission lines of AR Lacertae, estimated that with 
his (photographic) apparatus, only 8 per cent polarization would be the 
minimum detectable.’ It is possible that photoelectric methods may be 
able to reduce this value considerably —although the general problem of 
observing the equivalent width of a stellar absorption line photoelectrically 
has not yet been solved in practice. 

If a reliable and small lower limit to the actual polarization can be set 
from direct observation, then it would appear that certain detailed conclu- 
sions about the geometry of the star would be possible. From the symmetry 
of the magnetic field, it follows that a star viewed along the dipole axis 
ought to show no polarization. Intermediary directions will show polari- 
zation up to the equatorial maximum of 4 or 5 per cent. The direction of 
the polarization will always be such that the maximum equivalent width 
will be observed at right angles to the direction of the projected dipole 
axis, whether the star is viewed equatorially or not. 

If the not-unreasonable assumption is made that a star’s magnetic dipole 


and 
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and rotation axes are coincident, then we have here at least a theoretical 
basis for a determination of the direction of the projected rotation axis on 
the plane of the sky. 

* NotE: The work leading up to the writing of this paper was done at the Harvard 
College Observatory during the spring of 1950. 

! Babcock, H. W., Ap. J., 110, 126 (1949). 

2 Menzel, D. H., Jbid., 84, 462 (1936). 

3 Struve, O., [bid., 108, 155 (1948). 


HORMONAL INFLUENCE ON NUCLEAR SYNTHESIS. 1. 
ESTROGEN AND UTERINE GLAND NUCLEI 
By Max ALFERT AND Howarp A. BERN 
DEPARTMENT OF ZOOLOGY, UNIVERSITY OF CALIFORNIA, BERKELEY 


Communicated by Curt Stern, February 5, 1951 


In a recent publication! Salvatore has pointed out that injection of es- 
trogen into ovariectomized rats results in a doubling of nuclear volume 
in uterine gland cells. Progesterone does not produce this effect. Con- 
sequently the author concluded that his observations and other evidence 


“clearly show(s) the genic nature of rhythmical growth” and that ‘‘estro- 
gen acts on the genes by inducing their reproduction.” 

Inasmuch as these conclusions were based solely on karyometric analysis, 
it seemed desirable to test Salvatore’s hypothesis of genic duplication by 
comparing the amounts of some nuclear components in uterine gland nu- 
clei before and after doubling of nuclear volume. Recently developed 
photometric techniques*® * make it possible to measure amounts of protein 
as well as desoxyribonucleic acid (DNA) in single nuclei. Furthermore 
it has been shown repeatedly that the latter substance constitutes a reli- 
able measure of the state of ploidy of a nucleus—a definite amount of DNA 
being associated with a haploid set of chromosomes.‘ ~7 

For the purpose of this experiment several 4-month-old Long-Evans 5-1 
rats were ovariectomized; 6 weeks later, some of these were injected 
itramuscularly with 250 I. U. of estrone in sesame oil, others with 0.5 mg. 
of progesterone, daily for 4 days. Vaginal smears were taken at the time 
of sacrifice, 6 hours after the last injection, to determine the efficacy of 
castration and of hormone treatment. Uterine segments were fixed in 
acetic alcohol. Paraffin sections of the uterus from one castrate control, 
one estrogen-treated and one progesterone-treated animal were mounted 
together on slides and treated with the Feulgen reagent for determination 
of DNA and the Millon reagent for determination of tyrosine and trypto- 


{ 
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phane. The amount ot colored substance in single nuclei was determined 
from its specific light absorption measured with a sensitive phototube 
through a microscope. Techniques and methods of calculation are dis- 
cussed in a previous publication.” The results of the DNA measure- 
ments are given in table 1. 

The volume determinations incidental to these measurements confirm 
the observation of Salvatore, namely, that estrogen treatment caused a 
doubling of nuclear volumes in the uterine gland cells. Progesterone 
alone caused no such increase. The effect of simultaneous administration 
of estrogen and progesterone is still to be determined. The basic nuclear 
volume in the present material is slightly different from that reported by 
Salvatore (77.1 uv’ instead of 88.1 u*). This may be due to a difference in 
fixation and /or to differences in the strain of rats used and is not considered 
to be important. The Feulgen measurements clearly indicate that the 
nuclei in all three animals contain essentially the same amount of DNA. 


TABLE 1 
Amounts oF DNA (ARBITRARY UNITS) IN INDIVIDUAL NUCLEI OF UTERINE GLANDS 
OF THE RAT 
NUMBER 
OF NUCLEAR 
VOLUME AMOUN' 
ANIMAL MEASURED IN p® OF DNA 
Castrate control 20 .1 2.41 0.09 
Estrogen-treated 20 7.5 1.6 2.36 + 0.08 
Progesterone-treated 12 81.6 = 1.5 2.15 + 0.08 
Tetraploid. liver nuclei of 
a normal animal 10 5.04 * O15 
* Nuclear volume omitted, since there is often no consistent relationship among values 
obtained from nuclei of different tissues.* 


Furthermore, this amount is that of a typical diploid nucleus, as is demon- 
strated by comparison with values obtained from tetraploid liver nuclei of 
a control rat measured on the same slide. Thus, injection of estrogen has 
not resulted in doubling of DNA content in uterine gland nuclei, and poly- 
ploidy, as would be indicated by an increased amount of DNA, evidently 
has not been induced by the hormone. 

Doubling of nuclear volume, however, is a feature which normally ac- 
companies polyploidy within a tissue. In mouse liver, for example, three 
size classes and three corresponding classes of DNA content in a ratio of 
1:2:4 occur.4®7 It has been argued® that the increase in nuclear volume 
in such a case is probably not due to the increase of DNA alone but rather to 
increased protein content, inasmuch as such nuclei contain anywhere from 
8 to 15 times more protein than DNA." This argument is supported by 
observations on a number of cases, where changes in nuclear volume are 
not accompanied by changes in DNA content but by parallel changes in 


204 BIOCHEMISTRY: ALFERT AND BERN Proc. N. A. S. 


total protein content.* * |! Table 2 indicates that this is also true in the 
present case. 

The extinctions due to the Millon reaction are nearly the same in sec- 
tions of equal thickness of small and large nuclei; this indicates that the 
protein concentrations in nuclei of castrate control and estrogen-treated 
animals are very similar. The total amounts of protein in these nuclei 
are accordingly correlated with the nuclear volumes, and it appears that 
estrogen treatment has induced an approximate doubling of the protein 
content in uterine gland nuclei. 

This particular case of nuclear protein synthesis independent of DNA 
synthesis is similar to that reported by Schrader and Leuchtenberger.’ 
In the testis of the hemipteran insect Arvelius albopunctatus, spermatocyte 
nuclei in different lobes of the testis stand in a volume relation of 1:2:8 
but contain the same amount of DNA. The same concentration of pro- 
tein was found in the medium- and large-sized nuclei. The authors then 
draw attention to the amazing parallelism between nuclear size and pro- 
tein content and indicate their opinion that a rhythmic growth of extra- 
chromosomal proteins is involved. The functional significance of this 

TABLE 2 
PROTEIN CONCENTRATIONS IN UTERINE GLAND NUCLEI AS INDICATED BY THE EXTINC- 
TIONS (OpTICAL DENsITIES) DUE TO THE MILLON REACTION OveR 4 w Pato LENGTH 


EBLANK EcoRRECTED 
ANIMAL (12 Measurements) (10 M&asSuUREMENTS) (Test — BLANK) 


Castrate control 0.157 = 0.010 0.037 + 0.002 0.120 = 0.010 
Estrogen-treated 0 147 * 0.006 0.042 + 0.008 0.105 + 0.007 


phenomenon in Arvelius is obscure, since both types of spermatocytes 
eventually give rise to apparently normal sperm cells. In the estrogen- 
treated rat, on the other hand, the phenomenon of rhythmic protein syn- 
thesis may be related to the function of uterine gland cells which are pre- 
pared for secretory activity by estrogen. 

Protein measurements as employed herein do not permit distinction 
between chromosomal and extra-chromosomal proteins of the nucleus, 
but further attempts to specify the chemical nature of these proteins are 
being undertaken in order to obtain evidence bearing on the concepts of 
chromosome structure proposed by Mirsky and Ris.'* Exact duplication 
of one or all nuclear protein components, if confirmed by more refined 
methods, may indicate that a multiplication of some basic units is in- 
volved. Our investigations are being extended to include the study of 
other instances of hormonal alteration of nuclear and cytoplasmic constit- 
uents with special reference to the anabolic activities of the nucleus. 

Summary.-—Our work confirms that injection of estrogen into ovariec- 
tomized rats results in a doubling of nuclear volume in uterine gland cells. 
Photometric analysis shows that the increase in nuclear volume is paralleled 


\ 
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by an increase in protein content, while the amount of desoxyribonucleic 
acid remains unchanged. 

Acknowledgments.—This work was supported by research grants from 
the University of California Board of Research. We are greatly indebted 
to Professor A. W. Pollister, Dept. of Zoology, Columbia University, for 
allowing the senior author use of his laboratory facilities to conduct the 
measurements described herein. 
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THE STRUCTURE OF PROTEINS: TWO HYDROGEN-BONDED 
HELICAL CONFIGURATIONS OF THE POLYPEPTIDE CHAIN 


By Linus PAULING, ROBERT B. CorREY, AND H. R. BRANSON* 


SATES AND CRELLIN LABORATORIES OF CHEMISTRY, 
CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIAT 


Communicated February 28, 1951 


During the past fifteen years we have been attacking the problem of the 
structure of proteins in several ways. One of these ways is the complete 
and accurate determination of the erystal structure of amino acids, pep- 
tides, and other simple substances related to proteins, in order that infor- 
mation about interatomic distances, bond angles, and other configurational 
parameters might be obtained that would permit the reliable prediction of 
reasonable configurations for the polypeptide chain. We have now used 
this information to construct two reasonable hydrogen-bonded helical con- 
figurations for the polypeptide chain; we think that it is likely that these 
configurations constitute an important part of the structure of both fibrous 
and globular proteins, as well as of synthetic polypeptides. A letter an- 
nouncing their discovery was published last year.! 

The problem that we have set ourselves is that of finding all hydrogen- 
bonded structures for a single polypeptide chain, in which the residues are 


| 
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equivalent (except for the differences in the side chain R). An amino acid 
residue (other than glycine) has no symmetry elements. The general oper- 
ation of conversion of one residue of a single chain into a second residue 
equivalent to the first is accordingly a rotation about an axis accompanied 
by translation along the axis. Hence the only configurations for a chain 
compatible with our postulate of equivalence of the residues are helical 
configurations. For rotational angle 180° the helical configurations may 
degenerate to a simple chain with all of the principal atoms, C, C’ (the 
carbonyl] carbon), N, and O, in the same plane. 

We assume that, because of the resonance of the double bond between 
the carbon-oxygen and carbon-nitrogen positions, the configuration of each 

H Cc 


residue N—C is planar. 
No 


This structural feature has been 
verified for each of the amides that 
we have studied. Moreover, the 
rescnance theory is now so well 
grounded and its experimental sub- 
stantiation so extensive that there 
can be no doubt whatever about its 
application to the amide group. 
The observed C—N distance, 1.32 
A, corresponds to nearly 50 per cent 
double-bond character, and we may 
conclude that rotation by as much 
as 10° from the planar configuration 
would result in instability by about 
1 keal. mole~'. The interatomic 
distances and bond angles within 
the residue are assumed to have the 
values shown in figure 1. These 
values have been formulated? by 
consideration of the experimental 
values found in the crystal structure 
studies of pL-alanine,* L-threonine,* 
Dimensions of the polypeptide chain. N-acetylglycine®, and 8-glycylgly- 
cine’ that have been made in our 

Laboratories. It is further assumed that each nitrogen atom forms a hy- 
drogen bond with an oxygen atom of another residue, with the nitrogen- 
oxygen distance equal to 2.72 A, and that the vector from the nitrogen atom 
to the hydrogen-bonded oxygen atom lies not more than 30° from the N—-H 
direction. The energy of an N--H - - - O--C hydrogen bond is of the order 
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The helix with 5.1 residues per turn. 


x 
x 
2 
~ 


PAULING, 


CHEMISTRY: 


FIGURE 2 


Vor. 37, 1951 


The helix with 3.7 residues per turn. 
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of § keal. mole ', and such great instability would result from the failure 
to form these bonds that we may be confident of their presence. The 
N--H - - « O distance cannot be expected to be exactly 2.72 A, but might 
deviate somewhat from this value. 

Solution of this problem shows that there are five and only five configura- 
tions for the chain that satisfy the conditions other than that of direction of 
the hydrogen bond relative to the N--H direction. These correspond to 
the values 165°, 120°, 108°, 97.2° and 70.1° for the rotational angle. In 


the first, third, and fifth of these structures the eo group is negatively 


and the Ds ~H group positively directed along the helical axis, taken as 


the direction corresponding to the sequence—CHR—CO—NH--CHR 
of atoms in the peptide chain, and in the other two their directions are 
reversed, The first three of the structures are unsatisfactory, in that the 


FIGURE 4 
FIGURE 5 
Plan of the 3.7-residue 
helix. Plan of the 5.1-residue helix. 


N-—H group does not extend in the direction of the oxygen atom at 2.72 
A; the fourth and fifth are satisfactory, the angle between the N—H vee- 
tor and N—O vector being about 10° and 25° for these two structures 
respectively. The fourth structure has 3.69 amino acid residues per turn 
in the helix, and the fifth structure has 5.13 residues per turn. In the 
fourth structure each amide group is hydrogen-bonded to the third amide 
group beyond it along the helix, and in the fifth structure each is bonded to 
the fifth amide group beyond it; we shall call these structures either the 
3.7-residue structure and the 5.l-residue structure, respectively, or the 
third-amide hydrogen-bonded structure and the fifth-amide hydrogen- 
bonded structure. 
Drawings of the two structures are shown in figures 2, 3, 4, and 5. 
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For glycine both the 3.7-residue helix and the 5.1-residue helix could 
occur with either a positive or a negative rotational translation; that is, as 
either a positive or a negative helix, relative to the positive direction of the 
helical axis given by the sequence of atoms in the peptide chain. For 
other amino acids with the L configuration, however, the positive helix and 
the negative helix would differ in the position of the side chains, and it 
might well be expected that in each case one sense of the helix would be 
more stable than the other. An arbitrary assignment of the R groups has 
been made in the figures. 

The translation along the helical axis in the 3.7-residue helix is 1.47 A, 
and that in the 5.1-residue helix is 0.99 A. The values for one complete 
turn are 5.44 A and 5.03 A, respectively. These values are calculated for 
the hydrogen-bond distance 2.72 A; they would have to be increased by a 
few per cent, in case that a larger hydrogen-bond distance (2.80 A, say) 
were present. 

The stability of our helical structures in a non-crystalline phase depends 
solely on interactions between adjacent residues, and does not require that 
the number of residues per turn be a ratio of small integers. The value 
3.69 residues per turn, for the third-amide hydrogen-bonded helix, is most 
closely approximated by 48 residues in thirteen turns (3.693 residues per 
turn), and the value 5.13 for the other helix is most closely approximated 
by 41 residues in eight turns. It is to be expected that the number of resi- 
dues per turn would be affected somewhat by change in the hydrogen-bond 
distance, and also that the interaction of helical molecules with neighboring 
similar molecules in a crystal would cause small torques in the helixes, de- 
forming them slightly into configurations with a rational number of residues 
per turn. For the third-amide hydrogen-bonded helix the simplest struc- 
tures of this sort that we would predict are the 1|1-residue, 3-turn helix 
(3.67 residues per turn), the 15-residue, 4-turn helix (3.75), and the 1S-resi- 
due, 5-turn helix (3.60). We have'found some evidence indicating that 
the first and third of these slight variants of this helix exist in crystalline 
poly peptides. 

These helical structures have not previously been deseribed. In addi- 
tion to the extended polypeptide chain configuration, which for nearly 
thirty years has been assumed to be present in stretched hair and other 
proteins with the 6-keratin structure, configurations for the polypeptide 
chain have beeti proposed by Astbury and Bell,’ and especially by Huggins* 
and by Bragg, Kendrew, and Perutz.? Huggins discussed a number of struc- 
tures involving intramolecular hydrogen bonds, and Bragg, Kendrew, and 
Perutz extended the discussion to include additional structures, and 1n- 
vestigated the compatibility of the structures with x-ray diffraction data 
for hemoglobin and myoglobin. None of these authors proposed either 
our 3.7-residue helix or our 5.1-residue helix. On the other hand, we would 
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eliminate, by our basic postulates, ali of the structures proposed by them. 
The reason for the difference in results obtained by other investigators and 
by us through essentially similar arguments is that both Bragg and his 
collaborators and Huggins discussed in detail only helical structures with 
an integral number of residues per turn, and moreover assumed only a 
rough approximation to the requirements about interatomic distances, 
bond angles, and planarity of the conjugated amide group, as given by our 
investigations of simpler substances. We contend that these stereochemi- 
cal features must be very closely retained in stable configurations of poly- 
peptide chains in proteins, and that there is no special stability associated 
with an integral number of residues per turn in the helical molecule. Bragg, 
Kendrew, and Perutz have described a structure topologically similar to our 
3.7-residue helix as a hydrogen-bonded helix with 4 residues per turn. In 
their thorough comparison of their models with Patterson projections for 
hemoglobin and myoglobin they eliminated this structure, and drew the 
cautious conclusion that the evidence favors the non-helical 3-residue 
folded a-keratin configuration of Astbury and Bell, in which only one-third 
of the carbonyl and amino groups are involved in intramolecular hydrogen- 
bond formation, 

It is our opinion that the structure of a-keratin, a-myosin, and similar 
fibrous proteins is closely represented by our 3.7-residue helix, and that this 
helix also constitutes an important structural feature in hemoglobin, myo- 
globin, and other globular proteins, as well as of synthetic polypeptides. 
We think that the 5.1-residue helix may be represented in nature by super- 
contracted keratin and supercontracted myosin. The evidence leading us 
to these conclusions will be presented in later papers. 

Our work has been aided by grants from The Rockefeller Foundation, 
The National Foundation for Infantile Paralysis, and The U. S. Public 
Health Service. Many calculations were carried out by Dr. S. Wein- 
baum. 

Summary.—Two hydrogen-bonded helical structures for a polypeptide 
chain have been found in which the residues are stereochemically equiva- 
lent, the interatomic distances and bond angles have values found in amino 
acids, peptides, and other simple substances related to proteins, and the 
conjugated amide system is planar. In one structure, with 3.7 residues per 
turn, each carbonyl and imino group is attached by a hydrogen bond to the 
complementary group in the third anude group removed from it in the 
polypeptide chain, and in the other structure, with 5.1 residues per turn, 
each is bonded to the fifth amide group. 

* Present address, Howard University, Washington, D.C. 

+ Contribution No, 1538. 
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CONCERNING NON-CONTINUABLE, TRANSCENDENTALLY 
TRANSCENDENTAL POWER SERIES 


By F. BAGEMIHL 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ROCHESTER 
Communicated by J. L. Walsh, February 23, 1951 


The main purpose of this note is to show that power series of the kind 
described in the title can be obtained from a given power series by simply 
multiplying certain of its coefficients by — 1. 


Consider the class K of power series of the form }°a,2” whose circle of 


y=0 
convergence is the unit circle. There are ¢ elements in K (where ¢ denotes 
the power of the continuum). Let @ be the class of those series in K which 
can be continued beyond the unit circle, and let @ be the class of those 
series in K which satisfy an algebraic differential equation. Denote by 
e’, @’, the respective complements of ©, @, with respect to ®. 

There are the following sufficient conditions for a series in K to belong 
to @’, @’, respectively: 

(A)! Let {d,} (vy = 1, 2, 3, ...) be an increasing sequence of non-negative 
integers such that X,/v-> © asv-> », If Yaz’ belongs to ®, then it also 


vel 
helongs to @’. 
(B)? Let jdA,} (v = 1, 2, 3, ...) be a sequence of non-negative integers such 


that > vd, for every v. If belongs to ®, then it also belongs to @’. 


v=) 
The series }0s", which represents (1 — z)~! for | 
(1.e., to both @ and @). The series > 5,2", which represents the mero- 
morphic function ['(z + 1) for < 1, belongs to and’ to @’.  Accord- 


< 1, belongs to C@ 


ing to (A), $22” belongs to @’, and it is known‘ that this series belongs to 


(t. Finally, it follows from (A) and (B) that 3s” belongs to e@’@’. Thus, 


| 
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none of the classes C@, C@’, C’@, C’’ is empty. In fact, each of these 
classes contains ¢ elements; for if an arbitrary constant is added to any one 
of the four series just mentioned, the resulting series belongs to the same 
class as the original. 
Let }oa,2" belong to K. Call a sequence {e,} such that 6, = +1 (v = 
v=0 
0, 1, 2, ...) a “sign sequence.’’ Then® there is a sign sequence {e,} such 


that >¢,a,2" belongs to @’; indeed, there are ¢ such sign sequences. It 


is also known’ that there are infinitely many sign sequences {6,{ such that 


belongs to @’.5 We now prove the following 

Turorem. Let = f(z) belong to K. Then thereurec sign sequences 
v=O0 


fe,} such that >> ¢,a,2" belongs to a’. 
ve=O0 

Proof:* Since our series, by assumption, belongs to K, there is an in- 
creasing sequence of natural numbers {v,} (« = 1, 2, 3, ...) such that, for 
every kK, > a, # O, and such that lim | a,, | = lasx— ©. 
It follows then that © as According to (A) and (B), 
= g(s) belongs to @’@’; and the series a-g,,(z) = 
a=] 
where a # 0 is an arbitrary constant, {y,} is any infinite subsequence of 
and is any sign sequence, belongs a fortiori to C’@’. Set f(z)-g(z) 
= fo(z). Divide g(z) into an infinite sequence of power series /)(z), fo(z), 
...y fp(8), ..., each of which contains infinitely many terms, such that 
every term of g(z) belongs to precisely one f,(z) with p 2 1. Consider the 
set of all power series 


F(z) = fo(z) + 206,f,(z), 6, = #1. (1) 
y=1 


At most an enumerable number of these series can belong to @.° Hence, 
¢ of them must belong to @’, and these can be divided into ¢ pairs, because 
c=c+e. If Fi(s) and F,(z) are the members of any one of these pairs, 
then it is evident from (1) that F\(z) — Fi(s) = 2¢,,(s) for suitable se- 
quences » and a. As we remarked before, 2g,,(z) belongs to @’, so that"? 
at least one of F(z), F:(s) belongs to @’; and this completes the proof. 
The theorem does not remain valid if @’@’ is replaced by any one of the 
other three classes, because, according to (A) and (B), )Oe,2”, e.g., belongs 
to @’a@’ for every sign sequence {e,}, due to the presence of large gaps (i.e., 
consecutive terms whose coefficients are zero). There are, moreover, ¢ 
series in &, none of which has gaps, and yet each of which belongs to @' a’ for 
every sign sequence. For let {e,{ be an arbitrary sign sequence. Consider 


bed 
i 
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a 
264,2,, where ay) = 1, a, = v for every “belonging to the sequence {y,}, 
. 2 
where = 1,4%41) =¥% + «(« = 1, 2,3, ...), whereasa, = for every 


other v (this series obviously belongs to K). Every coeflicient is different 
from zero, and is an algebraic number, being either an integer or the recipro- 
cal of an integer. Furthermore, there is clearly no constant ¢ > O such that 


| = exp (log for every v 2 2. Consequently,'' 
v=@0 


belongs to @’. If v is not a term of the sequence |y,{, then €,a, ae 
The sequence }e,,a,,}, however, is unbounded, and — as 


©. Therefore’? cannot be bounded in any sector of the unit 
v=0 


circle, which means that this series belongs to @’.. Thus, }0¢,a,2" has no 
v=0 


gaps and yet belongs to @’@’ for every sign sequence {e,}. We can easily 
obtain ¢ such series, and, in fact, having rational coefficients, from Ya,2’ by 
v=0 
replacing a, = 1/v", for every » # 0 which does not belong to {»,}, by any 
2 
one of the numbers 1/77 (7 = v, vy + + 2,...). 


' Fabry, E., ‘Sur les points singuliers d’une fonction donnée par son développement 
en série et l’'impossibilité du prolongement analytique dans des cas trés généraux,” Ann. 
Sci. Ecole Norm. Sup., Ser. 3, 13, 367-399 (1896). 

2 Gronwall, H., “Sur les fonetions qui ne satisfont a aucune équation différentielle 
algébrique,”’ Ofversigt Kongl. Vetensk.-Akad. Forhand., Stockholm, 55, 387-395 (1898). 

8 Holder, O., ‘Ueber die Eigenschaft der Gammafunction keiner algebraischen Dif- 
ferentialgleichung zu gentigen,”’ Math. Ann., 28, 1-13 (1887). 

4‘ Jacobi, C. G. J., “Uber die Differentialgleichung, welcher die Reihen 1 * 2g + 
29? 4+ 2Wg + + + ete. Genitige leisten,” J. Reine Angew. 
Math., 36, 97-112 (1848). 

5 Hurwitz, A., and Polya, G., “Zwei Beweise cines von Herm Fatou vermuteten 
Satzes,"’ Acta Math., 40, 179-183 (1916). 

This follows from Hurwitz's argument, /hid., 182.183 

7 Ostrowski, A., “Uber Dirichletsche Reihen und algebraische Differentialgleichun- 
gen,”’ Math. Z., 8, 241-298 (1920). 

’ The existence of a sign sequence of this sort was recently rediscovered by Pélya, 
G., “Remarks on Power Series,’’ Acta Sct. Math. Sceged, 12, Leopoldo Fejér et Frederico 
Riesz LX X annos natis dedicatus, Pars B, 199-203 (1950). 

®» The proof makes use of ideas of Hurwitz, /. ¢., 182-188, and Ostrowski, /. ¢., 271. 

10 Moore, E. H., ‘Concerning Transcendentally Transcendental Functions,” Math 
Ann., 48, 49-74 (1897). 

4 Popken, J., ‘Uber arithmetische Eigenschaften analytischer Funktionen,”’ Diss 
Groningen, 1935; Theorem 12. A statement of this Theorem 12 can also be found in the 
more accessible article: Popken, J.,and Mahler, K., ‘Ein neues Prinzip fiir Transzend- 
enzbeweise,’’ Koninklijke Akad. van Wetenschappen te Amsterdam, Proc. Sect. Sci., 38, 
864-871 (1935). 

2 Duffin, R. J., and Schaeffer, A. C., “Power Series with Bounded Coefficients,” 
Am. J. Math., 67, 141-154 (1945); p. 158, Theorem IT. 
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SIMPLE WAVES IN THE COMPLETE GENERAL PROBLEM OF 
PLASTICITY THEORY 


By HILpA GEIRINGER 
DEPARTMENT OF MATHEMATICS, WHEATON COLLEGE 


Communicated by J. von Neumann, February 14, 1951 


The Problem.—The “complete” problem of the plane perfectly plastic 
body is determined by five equations for five variables: the components 
0;, 9,, 7 Of the symmetric plane stress tensor, 2, and the components, v,;, v,, 
of the flow velocity vector, 7; the system of equations is supplemented by 
adequate boundary conditions. One of the five equations, the ‘‘yield 
condition” is a finite relation between ¢,, o,, 7, stating that the plastic state 
of stress is restricted to a two-dimensional variety. 

In “plane strain,” ¢, = > “= 0, and for all other variables likewise: 

0/0z = UV. Because of the approximate incompressibility of the three- 
dimensional plastic medium one assumes, €, + €, + € = 0; hence e, = 0 
implies €, + ¢, = 0, and vice versa. With 
Ov, Ov, 


ay), 


Ox 


the five equations for ‘plane strain’ are! 


f(o,, 7) = Flo, 02.) = 0 “yield condition” (2) 


€:/0, = &,/o, = y/2r. (3) 


For reasons of symmetry f(¢,, o,, 7) in (2) depends, in case of an isotropic 
medium, merely on the two principal stress components, o;, 02. From (3) 


Ov Ox / Oy Ox 2 i 


(3”) 


we derive 


’ 


For ‘characteristic’? coordinates a(x, y) = const., B(x, v) = const., with 
Ug denoting the components of in the a- and 3-directions, (3") 
may be replaced by 


(4) 


= Ov Ox 
Ox Ox Oy 
Ov, Oy 
Ox Oy 
Ove , Og 
— v, = 0, v, = 0 
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equations which do no more contain explicitly the stress components. * 
In the literature of plane strain a particular case is almost exclusively 
considered: (2) is assumed as 


(o, — + 4r? — 4K? = — o2)? — 4K? = (2’) 


where A is the yield stress in shear. We shall call the problem (1), (2’), 
(3) the ‘classical problem” ;* we shall, however, see that many results usu- 
ally proved for (2’) are actually valid for general plane strain. 

In the typical plane problem of a thin plate under the action of forces 
applied to its contour parallel to the middle plane of the plate, the assump- 
tion, «, = 0, cannot be made. If then this assumption is abandoned, (3) 
must likewise be abandoned because (3) leads to (3"). Following von 
Mises, * (3) is then replaced by 


which for (2) = (2‘) reduces to (3). Equations (1), (2), (5) form the gen- 
eral complete plane problem. The ‘“‘stress problem” (1), (2) has been con- 
sidered only recently, by Sokolovsky,® von Mises,® Geiringer,’ Sauer® and 
Hodge.’ In the present paper we shall study the complete problem (1), 
(2), (5) with respect to basic equations, characteristics and initial value 
problems. Finally, for ‘‘simple waves’ both stresses and velocities will be 
explicitly determined. 

General Stress Problem.—TYhe three equations (1), (2) may be reduced in 
several ways to two equations. Sauer* uses, to advantage, the.two vari- 
ables o,, rt. Hodge® uses (0; + o2) and (0, — Following von 
Mises® we use the variables, s and 6, where @ is the angle of the x-axis and 
the first principal stress direction, and 5 is a parameter by means of which 
(2) is represented as = = o(s). Denote by &, the directions 
of first and second principal stress; (1), (2) are then replaced by 

O: 


08 Os 
fo (s) On On 


06 
g(s) de (6) 


where ae’ - are directional derivatives and, with = - 
ds 


f(s) = — (02 — 91), g(s) = (02 — 1). (7) 


The equations (6) are ‘“‘reducible.’’'! The plane of these variables, s, 4, 
is called ‘‘stress graph’ by von Mises. In this plane there are fixed char- 
acteristics, *, ~, determined by 
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d6 
+ or +6 = F(s) + const. (8) 


The characteristics, C+ C~ in the x, y plane, the physical plane, depend on 
the solution (on the boundary conditions) of a problem. We find for 
them :’ 

dy f(s). tan? g= f 


gg 


where doz/do, follows directly from the yield condition F(a, 02) = 0; also 


d 
m = i tan (9 = ¢). (10) 
dx 


OF 
General Velocity Problem.—From (2) with 5 = F,, (4 = 1, 2) we find 


oF OF OF 
= F, cos? 0 + Fy sin? 6, = F, sin? 6+ F, cos? 6, = 
Oo, 0a, Or 


(F, — F.) sin 20. (11) 


Because of (5), + e, = implies - = F, + F, = 0, and tan? ¢ = 
Oo, Oa, 


do./do, = —F,/Fo = 1, ¢ = w/4, C+ | C-. The characteristics, C+, 
form an orthogonal net, 1.€., ¢ = r if and only af €, + €, = 0, which in turn, 


is a consequence of (3); the C*, C~ then coincide with the shear lines; the 

general form of (2) ts then p(o, — o2) = O, or (without isotropy re- 

striction), 7 = g(a, — o,) where p and q are arbitrary functions of one variable. 

Ife + & + & = 0 ts assumed, this ts the case of “plane strain,” €, = 0. 
Using (5) as well as (9), (10) we obtain our basic velocity equations 


ov,  dv,\ //dv, cos2¢ //dv, dv, 
ov Ov Ov Ox sin 26 Ov Ox Oy Ox 
Ov, /OVv, 
— cot 26, = tan (6 + tan (@ — (12) 
Ox / Oy 
Any of the eqs. (12) follows from the two others. If the stress problem 
has been solved, i.e., if 6 = 0(x, vy), ¢ = ¢(x, ¥) is known, the eqs. (12) form 
a linear system; for g = m/4 they reduce to the particular eqs. (3’), (3”). 
It is seen that, as in the ‘‘classical’’ case, the characteristics of (12) coincide 
with the characteristics, C*, C~ of the stress problem (1), (2)." 
Let us introduce the v,, v,-plane, the velocity plane (Fig. 1). If in (12), in 
a domain where the respective Jacobien is not zero, we interchange de- 


{ 
216 

é 
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pendent and independent variables we find two equations whose character- 
istics, C+, C- exhibit a remarkably simple relation to the C+, C~: 


1 
—cot (0 = (13) 
m 


(13’) 


The relations between the C*, C” in the physical plane, *, I~ in the 
stress plane, C+, C~ in the velocity plane can be used to obtain in a very 


Phy si cal plane Velocity plane Speed draph 


FIGURE 1 


simple geometric way approximate solutions of the following classical 
initial value problems: 

(1) Along a nowhere characteristic curve, K, in the physical plane, s, 6 and 
Uz, Uy are given. 

(2) vw, (or v,) ts given along a Ct, v, (or v,) along a C~. 

Find s, 6, v,, v, in the characteristic quadrangle determined in (1) by K, 
in (2) by Ct, C~. We omit details. 

The eqs. (13) have a meaning independent of ‘‘velocity plane’ and ‘‘ex- 
change of variables.” They express for the eqs. (12) the relations which 
always exist between the derivatives of the dependent variables in a char- 
acteristic direction. Denote by d/d/*, d/dl~ derivatives in these direc- 
tions and put 29 = a — a’,a = 0+ 9, a’ = 6 — g; (13) may be writ- 


ten: 
dv, /dv, dv, /dv, 
- = —cot a’ 14 


From (14) we now derive the eqs. which generalize our basic eqs. (4). 
Denote by 1, 2, 3, 4 directions which make with the x-axis the angles a, 


= 
{ 
| 
dv, 
vy 
r* 
L 
P. 
- 
- 
0 x vy cr 
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a + 9’ a’, a’ + 9’ respectively, and by v, (¢ = 1, 2, 3, 4) the components 
of vin these directions. The result 1s: 


dv, da d Vs 


5 


If ¢ = w/4 (plane strain), 3 = 2, 1 = 4 and (15) reduces to (4). So far 
(15) contains, however, four components of v, and two definite components 
are still to be chosen. We now leave the general problem, and apply our 
equations to the case of simple waves where we can make a natural choice 
of appropriate components, and solve the whole problem explicitly. 

Simple Waves..The stress problem of simple waves has been considered 
by Geiringer.’ A simple wave!! region, D, of the physical plane is one 
that is mapped onto one characteristic, say I’) of the stress graph; to each 
point, s, @ on corresponds a straight characteristic, C*, in D, carrying 
constant values of s and 6; at all points of D (see (8)), 6 + F(s) is con- 
stant. A simple wave may be determined by giving I)” in the stress graph, 
and in the physical plane a set of straight lines, as the C*. It is then simple 
to find to each particular C*, or rather to its angle a, the corresponding 
constant solution, s, 6, of (6). 

We want now to determine the velocity distribution in D. Since a is con- 
stant along a Ct, eq. (15') with » = v,, and d/* = dr gives immediately 
(asin ‘‘plane strain’), 

dv, 
(16) 
dr 


Along a rectilinear characteristic C* the component v, of v in direction of C* 
remains unchanged. Next, from (15?) with = v, we derive 


dv, dv 
cos (a a’) — 


dl dl 


sin (a — a’) = [v, sin (a — a’) + 


da 
v, COS (a — a’)|}—. (17) 
dl 
Since the “cross characteristies,’’ C~, are not orthogonal to the C* we 
introduce as coordinate system the rectilinear C+ and their orthogonal 
trajectories. Denoting by T one particular trajectory we use as coordinates 
(a, r) of a point, P;, in D, the angle, a, which the C,* through P; makes 
with the x-axis, and the directed distance, 7 = OP, where O is the point of 
intersection of C;* and 7. Then, if ds, denotes the line element of the 
orthogonal trajectory in P; and R the radius of curvature of 7 in O we 


find that ds,,/da = r — R, and from (17) we obtain 


i 

i 

, 
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ou, ov 
cot 29 + (r — R) (18) 


which for ¢ = 7,4 reduces to v, + > = 0 (see (4)). Since R, a, ¢ do not 
a 


depend on r we find from (18) 


(r — R) cot 29 (19) 


an eq. of second order for v, alone. [t is, however, simpler to write (18) as: 


(R — r) + tan 2 +v, = ~ tan29= (a). (20 
a (Va) wla ) 


Ou Ov 
Along a C* both v, and w(a) = (R — r) > *+ tan 29 5 “+ v, are constant. 
r a 


It is now easy to solve the Cauchy problem, 1.e., to find the velocity dis- 
tribution in a simple wave region, D, if the velocity, v, ts known along a non- 
characteristic curve, K, across D. Assume along A, v, = h(a), and v, = 
k(a) as given. From k(a@) and from the expression of tan 2g in terms of 
a (which is known if a particular yield condition and a particular I'y~ have 
been chosen) we find the w(a) of (20) and the problem reduces to the solu- 
tion of the linear first order equation (20). Its characteristics are 


du, 


= = » 
R-r tan 29 — (21) 
Each of these ordinary equations contains only two variables, and is linear. 
The first dr da = (R — r) cot 29 is the eq. of the cross characteristics, C 
If fila, r) = a, fola, r) = are the solutions of (21), then G( fi, = O 
(where G is an arbitrary function) is the general integral of (20). To adapt 
the so found v,(7, a) to the initial condition, v, = A(a) along AK, we have to 
express the r along A in terms of a. Details of the procedure, as well as 
examples must be omitted. The ‘‘characteristic initial value problem” 
(see p. 217, (2)) can likewise be solved explicitly. 

For the general complete problem (1), (2), (5) we have thus obtained 
the following results: (a) the velocity equations (12), (14), (15); (db) 
approximate solutions of the classical initial value problems; (c) exact 
simple wave solutions. Solutions of initial value problems form the ma- 
terial for building up solutions of boundary value problems. 
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§ For more details see: (a) Geiringer, H., and Prager, W., Ergeb. exakt. Naturw., 13, 
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i] 
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Oa 
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CONTINUOUS GAMES 
By SAMUEL KARLIN* 
CALIFORNIA INSTITUTE OF TECHNOLOGY AND PRINCETON UNIVERSITY 


Communicated by J. von Neumann, February 28, 1951 


In this note a new method is presented for handling continuous games. 


One can deduce both new theoretical results and also qualitative statements 
as to the nature of optimal solutions for certain classes of kernels arising 
in continuous games. The essential idea is to vary the spaces of strategies 
and study the effect this has on the game. The methods employed involve 
partial orderings in abstract spaces. Throughout, we use the usual ter- 
minology of the Theory of Games. 

1. Let & and y denote Banach spaces consisting of measures and let ® 
and § denote generally spaces of functions such that ®* and $* (conjugate 
Banach space) consist again of measures. Let A denote a closed cone of 
elements in @, i.e., if x, ve A, then Ax + uve A ford, uw > 0, andifxe K 
and —x isin A, then x = 0. Similarly, P shall denote a closed cone in &. 
The dual cone A* in &* shall be defined as the set of all fin (* with the 
property that (x, f) > O for every x in A. [(v, f) is the value of the fune- 
tional f acting on the element x. Throughout this paper the inner product 
notation shall be used.| With respect to any cone one can define a partial 
ordering, e.g., x > Oif and only ifx «A. Similarly y > x shall mean y — 
x > 0. Also, let Q be a closed cone contained in A*. We assume that the 
cone A possesses an interior point «. A and B shall denote bounded linear 
operators mapping & into ®, and the conjugate mappings shall take R* = 
y into S. A; shall denote the section of the cone A*, where (u, f) = 1. 
Q,, is defined analogously with the section taken relative to Q. Finally, it 
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is assumed that in 8 a cone L is given and that P is contained in L*. Also, 
L has an interior point v to which is associated P,. 

Definition 1: Let Q,(A, B) be the set of real numbers \ for which there 
exists a fin P, such that 


Af > \Bf (Af — ABf isin 


Put A, = B) = sup A for in Q,. 

Definition 2: d, and ®, are defined in a similar manner with the above 
inequality reversed and sup replaced by inf. 

Definition 3: Let Q,(A, B) = [| for any gi, ..., € Q, there exists a 
fin P, with (g,, Af) > Agi, Bf)t = 1, ..., nm]. Clearly, f may depend 
upon g,. Put), = forA in Q,. Similarly, &, and A, are defined. 

Definition 4: It is furthermore assumed that P and Q are chosen large 
enough. More precisely, if (g, x) > O for every g in Q,, then x is in K. 
The set Q, is said to be sufficient if whenever x is of the form x = Af — 
\Bf for fin P, and (g, x) > O tor every ginQ,, thenxisin A. Ananalogous 
statement applies to P. 

Definition 5: 0,(A*, B*), 2,(A*, B*) ... and = 4,(A*, B*) ... are 
defined analogously in terms of the adjoint mappings and Q,,. 

The meaning of these quantities is explained in the following result: If 
B is assumed to be strictly positive in the sense that for any g ¢ Q, and f ¢ 
P, it implies (g, Bf) > 0, then one gets with sufficiency that 


sup inf = (g, Af) 
feP, (g, Bf) 


inf sup (& Af) _ 
g€Q, feP. (g, Bf) 


The assumption on B is not really essential but is made here to avoid un- 
necessary writing. These formulae indicate the connection with the funda- 
mental quantities occurring in the Theory of Games. The general relation- 
ship between these quantities is given in the following Theorem. We as- 
sume the requirement of sufficiency ts fulfilled. 

TuHeoreM \, > AT >A, > AL 

In the sequel, the statement that a value exists shall be synonymous with 
the equality of A = ),. It is to be understood that the value exists means 
here relative to the cones P and Q. 

Any one of the following conditions will insure At = ),. 

THroreM 2. (a) If the image of Af — X Bf for { in P, and any fixed 
constitute a weakly compact set relative to Q,,, then this implies that \, = Xz. 

(b) If the image of A*g — X\B*g with g in Q, 1s weakly compact relative to 
P,, then = 

Remark: A set / of elements is said to be weakly compact relative to a 


| | 


222 MATHEMATICS: S. KARLIN Proc. N. A. S. 


‘set G of functionals if under the weak topology on // induced by G, then 
I] constitutes a bicompact set. 

Many other slight variations of the conditions of Theorem 2 can be given 
which will imply the existence of a value. The advantage of this general 
formulation is the great flexibility in choosing the cones. It is clear that 
if PY and ¢ imply < AO < Thus to 
secure a value one wants to increase the cones P and Q sufficiently until one 
obtains A, = A*. On the other hand, if one has a value, then to obtain 
more refined qualitative results it is desirable to reduce the size of the cone 
to a minimum cone which still has a value. Furthermore, keeping Q, 
fixed and decreasing P, still maintains a value as long as P is a cone suf- 
ficient in the sense of Definition 4. Some applications of these ideas will 
be presented in Section 2. Also an equivalence theory for solution can be 
developed based on this principle. Two cones P“ and P™ are said to be 
equivalent relative to the operators A and B if \,(P") = \,(P*). 

Finally, in this section some results on perturbation of games and its 
applications are stated. 

THroreM 3. (a) If A, = A (uniformly) and B,, = B, then \" and 
The cones P and Q are held fixed. 

(b) If Ay 2 Az => ... ~* A (strongly) and A* — A* strongly with B, = B, 
then if A, maps & into Rand & ts the conjugate space to the space containing 
the image of A}, then A has avalue. The cones P and Q are held fixed. 

(c) If P! ¢ P? P (strongly) while also increases to Q and if 
with respect to each P?, Q% a value exists, then a value also exists relative to 
the cones P and Q. 

A few examples will illustrate this theory. 

(1) Let Af = {/M(x, y)f(x)dy, where M(x, y) is a measurable function 
of both variables and bounded on the unit square. Let B = (v, /)u, where 
u = v is the function identically 1. Finally, let P represent the cone con- 
sisting of all absolutely continuous distributions over the imterval (0, 1). 
With the aid of Theorem 3 it can be shown if Q = P, then a value exists. 
In this case no optimal solutions in general exist, however, the value may be 
approached within ¢ by suitable choice of strategy. 

(2) If M(x, y) is a lower semicontinuous function, then a value exists 
relative to all completely additive distributions on the interval (0, 1). The 
definitions of A and B are taken to be the same operators as in the previous 
example. This example was first achieved by I. Glicksberg. 

2. In this section we indicate how qualitative results may be obtained 
by selecting the cones appropriately. The definitions of A and B are given 
as in example (1) above. The kernel J/(x, y) is in addition assumed to be 
measurable in each variable for fixed choice of the other variable. We dis- 
cuss several examples. 

(1) Let M(x, vy) = O with V,, being continuous in both variables. 
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In other words, M(x, y) for each x.is convex in y. We choose for P the 
closure of the set of all distributions which have only a finite number of 
points of increase. The space ® is taken to be the set of all functions 
possessing continuous second derivatives. The cone A is chosen to be all 
functions for which w"(y) > 0, w'(0) > 0, w(0) > 0. The condition Af > 
\ Bf is equivalent with 

So! y)df(x) > 0 

So! M(x, O)df(x) > 0 

M(x, O)df(x) > X. 
The first inequality imposes no restriction on any distribution. This fol- 
lows from the hypothesis put upon .J/(x, vy). Thus it may be discarded and 
it can be concluded from the last two inequalities that an optimal solution 
distribution for player | exists using at most two jumps. 

(2) If AV,...,(x, y) 2 O with the derivative taken n times, then an opti- 
mal solution distribution for player | can be found with at most m points 
of increase, 

and N,(x, vy) > 0. 
t 0 
Again, let A be chosen as the set of all continuously differentiable functions 
with w’(v) > Oand w(0) > 0. Let Q be chosen to be all possible absolutely 
continuous distributions spread over the set TP of y where N(y, y) < 0 and 
a possible jump at the origin. The condition Af — ABf > 0 becomes 
N,(0O, 9) + So” Ni (x, + Noy, vif(y) > 

aN(0,y) -. 
The first inequality is automatically satisfied over the set of y where 
N(y, y) 2 0. Thus the set Q designated above becomes sufficient for this 
problem. Therefore, it can be shown that within ¢ an optimal solution for 
player 2 can be found in the cone Q. 

Other types of kernels of this kind with a discontinuity along the diagonal 
can be handled in such a manner. 

(4) If M(x, y) is analytic for each x in y, then a suitable ordering for A 
is to require that w'"’(0) > 0 for every n. Special results can be obtained 
in this circumstance. 

These examples serve to illustrate the method. In essence, the idea is 
to utilize the properties of the kernel (x, y) and to accordingly select a 
cone within which a solution can be found. 

3. The entire above procedure and method can be applied to the theory 
of linear programming. There the concept of product operators and prod- 
uct spaces are used. Moreover, many applications of the above theory can 
be found in the theory of Decision Functions and the problems of testing 
hypothesis. Details of the above result and additional applications will be 
published later. 


(3) Let M(x, y) = 


*O. N.R. Fellow. 
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THEORY OF TURNS AND SLIDES UPON A SURFACE 
By EpwWArpD KASNER AND JOHN De Cicco 


DEPARTMENTS OF MATHEMATICS, COLUMBIA UNIVERSITY, NEW YorRK, AND DE PAUL 
University, CHICAGO 


Communicated February 7, 1951 
1. Consider a surface } for which the linear element is 
ds? = K(x, y)dx* + 2F(x, y)dx dy + G(x, y)dy? (1) 


such that / > 0,G > 0, H > 0, F, H® = EG — F* > O are real analytic 
functions over a certain region of the parametric (x, y)-plane. Thus & is 
represented in a point-to-point fashion upon this region of the plane. 

A dilatation D, on & is a lineal element transformation whereby any lineal 
element (x, y, P) corresponds to a lineal element (X, Y, P), such that both 
are orthogonal to the same geodesic through the points of the elements, 
and the geodesic distance between the given points is the constant number s. 

The set of dilatations D, on » form a one-parameter group. The infini- 
tesimal lineal element transformation, which generates this one parameter 
group, is the infinitesimal dilatation D,,. 

The infinitesimal dilatation D,, is given by the equations 


H(E + 2Fp + Gp?) ix = (F + Gp)és, 
H(E + 2Fp + by = —(E + Fp)bs, 


1/ 
+ 2Fp + Gp*)”* bp | (; + p 2) (E + 2Fp + 
Ox oy 


Gp*) + + 2Fp + Gp’) (= +p °) és. (2) 


An infinitesimal dilatation upon 2 is an infinitesimal contact transfor- 
mation with the characteristic function 


W = 2 (E + 2Fp + Gp?) (3) 
The dilatations upon Y form a one-parameter group of contact transfor- 
mations whose path curves are the geodesics of 2. 

Two infinitesimal contact transformations, neither of which are extended 
point transformations, such that they possess the same transversality law and 
their commutator is an extended point transformation, which 1s not the identity, 
are two infinitesimal dilatations of two conformally equivalent surfaces. 

Two infinitesimal contact transformations possess the same linear in- 
volutorial law of transversals if and only if they are two infinitesimal dila- 
tations of two conformally equivalent surfaces. ' 

2. Aturn upon a surface is a lineal element transformation whereby 
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each lineal element is rotated about its point through a constant angle a. 

The turns 7’, form a one-parameter group generated by the infinitesimal 
turn 7;,, namely 

6x = 0, by = 0, Hip = + 2Fp + Gp*)éa. (4) 
The path series are the ©? stars, each consisting of the ©! lineal elements 
through a given point of 2. 

A slide S, upon a surface & is a lineal element transformation whereby, 
to each lineal element (x, y, ») of 2, there corresponds a lineal element 
(X, Y, P) of X, such that it is tangent to the geodesic of © determined by 
(x, y, p), and the geodesic distance between the points of these given lineal 
elements is the constant k. 

The slides S, form a one-parameter group generated by the infinitesimal 
slide namely 
(E + 2Fp + = dk, 

(E + 2Fp + Gp’)'“éy = pok, 
+ 2Fp + = [(FE, — 2EF, + EE,) + 
(GE, — 2EG, — 2FF, + 3FE,)p + (—3FG, + 2GE, — EG, + 
2FF,)p? + (—GG, + 2GF, — FG,)p* 6k. (5) 
The path series are the * geodesics. 

3. If U,F, U,F UF, denote the symbols of the infinitesimal turn, slide, 
dilatation, then their commutators are related as follows: 

(Ui, U2) F = (U2, Us)F = KU,F, (U3, Ui) F = U2F, (6) 


where K denotes the Gausstan curvature. 


By means of these formulae, the following result 1s obtained. 

The turns and slides upon a surface = form a three-parameter group, called 
the whirl group, if and only if X is of constant Gaussian curvature. 

If a surface is of variable Gaussian curvature, then the turns and 
slides upon ¥ do not generate a group of any finite number of parameters.” 


' Kasner, Differential Geometric Aspects of Dynamics, Princeton Colloquium, T. Am. 
Math. Soc., 1913, 1949; Lie, “Die infinitesimale Bertihrungstransformationen der Me- 
chanik,"’ Leipziger Berichte, 145-153 (1889); Kasner, “The Infinitesimal Contact Trans- 
formations of Mechanics,” Bull. Am. Math. Soc., 16, 408-412 (1910); De Cicco, “New 
Proofs of the Theorems of Kasner Concerning the Infinitesimal Contact Transformations 
of Mechanics,” J. Math. Phys., 26, 104-109 (1947). 

2 The geometry of turns and slides was initially studied by Kasner in 1911, followed 
by a series of papers by Kasner and De Cicco. Recently, contributions to this subject 
have been made by A. Narasinga Rao, K. Strubecker, and J. M. Feld. See Kasner, 
“The Group of Turns and Slides and the Geometry of Turbines, Am. J. Math., 33, 193- 
202 (1911); De Cicco, ‘The Geometry of Fields of Lineal Elements,’’ Trans. Am. Math. 
Soc., 47, 207-229 (1940). This article contains a bibliography of the subject. Feld, 
“4& Kinematic Characterization of Series of Lineal Elements in the Plane and of Their 
Differential Invariants Under the Group of Whirl-Similitudes and Some of Its Sub- 
groups,” Am. J. Math., 70, 129-138 (1948). 


MATHEMATICS: K. MENGER Proc. N. A. §S. 


PROBABILISTIC GEOMETRY 


By KARL MENGER 
ILLINOIS INSTITUTE OF TECHNOLOGY 
Communicated by Marston Morse, February 11, 1951 


We genera = the concept of a metric space by associating a cumulative 
distribution function, A,,, with every ordered pair (a, 5) of elements of a 
set S. The value A,,(x) may be interpreted as the probability that the 
distance froma tobbe <x. But S need not be a metric space in the ordi- 
nary sense. The distribution functions and the association of these func- 
tions with the pairs of elements of S, are all that is assumed. For instance, 
a probabilistic metric consisting of five elements is a 5-by-5 matrix of dis- 
tribution functions as an ordinary metric space consisting of five points 
isa 5-by-5 matrix of numbers. 

In analogy to the assumptions 1|’., 2’., 3’. about a metric space,* the non- 
decreasing functions 4,,, which are continuous to the right, are assumed to 
satisfy the following postulates. (Probabilistic interpretations are men- 
tioned parenthetically. ) 

Ly”. Maal) = 1 for every point a and every x > 0. (The distance from 
a point to itself is certainly <0.) 

1,”. A(x) = O for every two points and every x < 0. (Negative 
distances are impossible. ) 

2”, Aagn(x) = Aya(x) for every two points and every x. (The distance 
from a to b is as likely to be <x as is the distance from 6 to a.) 

The triangle inequality included in our original set of postulates! will be 
replaced by the following more stringent version due to Wald.” 

3”. [Ag Ane < Age(x) for every three points and every x. Here, 
A’ » A” denotes the convolution of A’ and A”, that is, the distribution of 
the sum of the random variables with the distributions A’ and A”. [For 
every x, it is less likely (or, at any rate, not likelier) that the sum of the 
distances from a to 6 and from } to c be < x, than that the distance from 
atocbe < x.] 

An ordinary metric space with the distance function d(a, 6) is a special 
case of a probabilistic metric, namely, the case that 


soifx < d@, d) 
= tx > de, 


Normal distributions are ruled out by Postulate 1,”. 

We call a and a’ certainly-indistingutishable if Aga(x) = 1 for each x > 0. 
Uniting all elements which are certainly indistinguishable from each other 
into identity sets we decompose the space into disjoint sets A, B,... . 
We may define Ayg(x) = A,(x) for any a belonging to A, and any 4 be- 
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longing to B. (The number,is independent of the choice of a and 6.) The. 
identity sets form a perfect analog of an ordinary metric space since they 
satisfy the condition 

1.”". If A # B, then there exists a positive x such that Ayg(x) > 1. 

Replacing, if necessary, points by identity sets we shall assume that no 
two points are certainly-indistinguishable unless they are identical. 

Let 6(a, b) denote the greatest lower bound of all numbers x such that 
Aw(x) = 0. If this certainly-lower-distance is positive, then a and b will 
be said to be lf 6) = we distinguish two 
cases: we call a and 6 barely-distinguishable if A,,(+0) = 0; and perhaps- 
indistinguishable if Aw(+0) > 0. Here A,»(+0) denotes the right-side 
limit of A,, at 0. 

The relation of being perhaps-indistinguishable is reflexive, symmetrical 
and (by virtue of 3”.) transitive. Again the space can be decomposed into 
disjoint lumps consisting of perhaps-indistinguishable elements. But if 
a and a’ belong to the same lump, and + is any element, the functions Ags 
and A,,, need not be equal. Even if the distances between lumps are de- 
fined by distribution functions rather than numbers, the definition of these 
functions must be based on averaging processes. 

Now let S be an ordinary metric space with the distance function d(a, }). 
We call a probabilistic metric monotonic on S if d(a, b) < d(a, ©) implies 
Mac(a) < Ag(x) for every a, 6, ¢, x. Clearly, in a monotonic metric, 
d(a, b) = d(a, c) implies Ag, = 

If the metric space is complete, convex and externally convex® (so that 
every two points are on at least one straight line), then it can be shown 
that either each lump contains only one point or that there is only one 
lump, namely, the entire space. More precisely, we obtain the following 
theorem. One of three cases is present: etther every two points of S are per- 
haps-indistinguishable; or every two points are barely-distinguishable; or 
every two points are certainly-distingutshable. 

Each of these three cases can occur on the straight line, even in conjunc- 
tion with the assumption of 

Homogeneity. = (x) for every b, x. 

In a homogeneous metric of the line it 1s sufficient to define Ag,(x) fora 
> Oandx > 0. 

With the following quasi-exponential distributions, every two points are 
perhaps-indistinguishable. 

Aoa(x) = 1 — (2a > 0,x > O). 
Here, ® must be a nowhere decreasing function for which 
< 2b(a) — P*(a). 
A solution satisfying the postulates is ®(a) = | — e 
Exponential distributions, = 1 — e 


a 


, are impossible on the 
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straight line (and in much more general spaces) even if the assumption of 
homogeneity is dropped. 
With the following rectangular distributions, every two points are 
barely-distinguishable. 
Aug(x) = x/g(a) for 0 < x < ¢(a), 


where g(a) + g(b) — g(a + b) > gla)y(b)/2¢(a + &). A solution satis- 
fying the postulates is g(a) = a* if 
S14+ 1/72 or 0< 0.7714... 
A particularly simple example is g(a) = Va. 
With the following rectangular distributions, every two points are cer- 
tainly-distinguishable. 
x — ¥(a) 


Aoa(x =e 
g(a) — 


for Y(a) < x < ¢(a), 


where 
[e(a) — ¥(a)] — 
2[e(a + b) — + 


The postulates are satisfied if we set g(a) = ak (0<k < 0.7714...) asin 
the preceding example, and ¥(a) = cy(a) for some constant ¢c < 1. 

The metrics of the last two examples contract the scale in the sense that 
the distance from 0 to a is certainly <¢g(a), for instance, < Va. We are 
reminded of the psychophysical contractions of the scale, especially of 
Cramer’s suggestion’ that the subjective value of an amount of money, a, 
does not exceed 4/a. The assumption that the distance from 0 to a be 
certainly < log a (that is, the analog of the well-known Bernoulli-Fechner 
contraction) contradicts the triangle inequality fora < 2%. 

An inhomogeneous metric, which does not contract the scale is defined 
by the following rectangular distributions 


y(a) + ¢g(b) — g(a + b) > 


Anca) = Aeo(x) = Ao, for b < ¢ and every x, 
Aog(x) = x/afor0 <x <a. 


In this metric, every two points are barely-distinguishable. 

For some purposes, none of the three main cases is desirable. For in- 
stance, none of them is compatible with the existence of a point a, which 
lies in a bounded threshold consisting of points that are perhaps-indistin- 
guishable from a, while the points outside are certainly (or at least barely) 
distinguishable from a. From our theorem it follows that, in order to 
cover this and similar situations, we have to weaken, if not to give up, the 
triangle inequality. 

It appears that in the transition from the rigid geometry in the large to a 


) 
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probabilistic geometry in the small, the geometric laws to which we are accus- 
tomed, not only cease to be certain but, in some cases, become unlikely and even 
im possible. 

There seem to exist two or three main types of applications of a prob- 
abilistic geometry. 

1. To the psychophysical continuum. The observer starts with an 
objective metric space (or what, in a first approximation, he considers as a 
metric space): the skin of the blindfolded observed person, the sets of 
weights read on a scale, temperatures, and the like. Upon these objective 
rigid metrics, the observed person superimposes his subjective probabilistic 
metric. 

2. To the physical macrocosm. Our theory provides us with a scheme 
for dealing with bundles of observations of spatial quantities, weights, tem- 
peratures, etc. Such a scheme should be useful where the enormously 
simplifying method of replacing these bundles by single numbers would 
prove to be an oversimplification. Underlying the (more or less arbi- 
trary) averaging and selection procedures of the simplifying method is the 
idea that various phases of the macrocosm have the structure of a numeri- 
cal continuum (for the elements of which we have to discover the ‘real’ 
values and distances). But the foundation of this idea seems to lie only 
in scientific usage. 

3. Perhaps to the physical microcosm. In this application, the funda- 


mental difficulty seems to lie in the lack of an objective metric comparable 
to that of the observer of another person's psychophysical continuum. 
The postulate that each world point can be characterized by four numeri- 
cal coordinates is not more plausible than the assumption that the dis- 
tance between two points be a number. Consequently, the theory would 
have to develop criteria for identifying the elements of a set, along with 
the distance distributions between the elements. 


* See my note “Probabilistic Theories of Relations”’ in jthe preceding issuc of these 
PROCEEDINGS, 37, 178, 1951, where the assumptions 1.’, 2.’, 3.‘ are quoted on p. 179. 

! Proc. NaTL. ACAD. Sct., 28, 535 (1942). 

2 [hid., 29, 196 (1943). 

3 “Untersuchungen ueber allgemeine Metrik,”’ Math. Ann., 100, 75-163 (1928). 

4 Cf. the author’s paper in Z. Nationaloekonomie, 5, 465 (1934). 
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NUMERICAL RESULTS FOR PRESSURE-VELOCITY CORRELA- 
TIONS IN HOMOGENEOUS ISOTROPIC TURBULENCE 


By D. NELSON LIMBER 
YERKES OBSERVATORY, WILLIAMS Bay, WISCONSIN, UNIVERSITY OF CHICAGO 
Communicated by J. von Neumann, February 8, 1951 
Introduction.—In a recent paper Batchelor! has considered fluctuations 
i pressure in homogeneous isotropic turbulence by studying the correla- 
tion, pp’, of the pressures p and p’ at two points separated by a distance 


r. He has shown that this correlation, which in homogeneous isotropic 
turbulence is a function only of 7, can be expressed in terms of the scalar, 


Q, defining the fundamental correlation term u,v,’ if the assumption is made 
that the fourth order correlation u;4ju;'U»’ is related to the second order cor- 


relation uu,’ as if uy and were jointly normally distributed. 
Batchelor’s final result is that 


It is the object of this note to show that with the same assumptions as 


Batchelor has made, the correlation Puy’ u,’ can also be expressed in terms 
of Q. 
Theoretical Development.— Multiplying the Stokes-Navier equation 
Ou, ou, 1 op 
‘4+ = *u 
ot Ox, p OX, 
by u,'u,,' and averaging we obtain 


OP im , 


Uy Um" vv = dé, (uuu su,’ ttn’) + 


where £; = x,;/ — x, and where P,,, denotes the correlation tensor Pu,'u,,’. 
Under conditions of isotropy Pm must be of the form 


Pin PritiEn + Prbim, (3) 


where P; and P2 are two even functions of the distance r separating the 
points considered. 

On the assumption that the fourth order correlations are related to the 
second order correlations in the manner of a joint normal distribution, 


= QuQym + + (0)Om(0) (4) 


where it is known that? 


VoL. 37, 1951 PHYSICS: D. N. LIMBER 231 


Oy = ua! = © 10" + 20)by. (5) 


Since the quantity on the left-hand side of equation (2) is solenoidal in 7 
and symmetrical in / and m, the same must be true of the terms on the 
right. The term 


X int = dE, + QimQyi) (6) 


must therefore be solenoidal in 1 and symmetrical in / and m, in other 
words it must be derivable in terms of a single scalar XY in the manner® 


Nime = curl X(&; emia Ee + Em = 


(rX" + 3X) + Emus) + (7) 


The condition that equation (6) be of the form of equation (7) suffices to 
determine P; and P, in terms of Q. We find 


rP, + 3P; = —4 dy (8) 


+ ay (9) 


r 


24 1 14 
Py = 5 (Q')? dy + 5r8 (Q’)*y* dy. (10) 


The derivation of these results follows the development contained in a re- 
cent paper by Chandrasekhar on turbulence in magneto hydronamics.* 
Equation (8) is of special interest since by equation (3) 


1 —- 
Pilr) = r?P, + 3P2. (11) 


— Palr) = —4 So’ (Q')*y dy. (12) 


Results.It is of special interest to see what form equation (12) takes 
when the Kolmogorof Region is considered. We then have‘ 


Qu? € aS) = (ve) (")" (13) 


Thus 

and 

1 _ 

(14) 
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Subsututing the latter into equation (12) we obtain 


or by equation (13) 
(16) 


For the correlation pp’ Batchelor! has obtained the following relation 
under the same conditions: 


l l r s/s 


‘Thus we have the simple relation 


The longitudinal and transverse correlations of P,,, are, respectively, 
p(0, O)uz(x, O)us(x, 0) = Py(r)r? + P2(r). (19) 
and 
(0, 0)uy(x, O)u,(x, 0) = Py(r). (20) 


Now Proudman® has tabulated Q and Q’ for Reynolds numbers 12.5 and 
©, His tabulations are in terms of a dimensionless variable, y, given by 


r= y= oy, (21) 


where a is the R of Chandrasekhar,’ v is the kinematic coefficient of vis- 
cosity, and ¢ is the time elapsed since the beginning of the decay. In 
terms of the function f — g tabulated by Proudman, expressions (8), (10), 
and (19) take the form 


*P\(y) + 3P. 1 


= +> +5, f (f — g)txtde, (23) 
u 5 Sy x Jo 


ya*Pi(y) + Poly) (24) 


Using Proudman’s tables for infinite Reynolds number, the integrals have 
been evaluated and the results are presented in table 1. From these values, 


uf 
and 
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the correlations of expressions (19) aad (20) have been plotted in figures 
1 and 2. 


ORDINATE 
+B) 
u 


i 
2 3 
FIGURE 2 


5.0 


~P; 108 57. Q +0) 64 


+ 3P2) X 10? - 72 —2.4 - 0.68 
u 


u 


1 
= + P2) X10? 17 +0.10 +1.07 +1.18 +1.04 +0.84 
u 


It should be possible to check these predicted correlations experimentally 
and thus obtain a better understanding concerning the underlying theory. 

In conclusion I should like to express my gratitude to Professor Chan- 
drasekhar who not only suggested the investigation but also has given gen- 
erously of his time and patience in seeing it to its completion. 

! Batchelor, G. K. (in press). I am indebted to Dr. Batchelor and to Professor 
Chandrasekhar for having made a manuscript copy of this paper available to me 

2 Robertson, H. P., Proc. Camb. Phil. Soc., 36, 209 (1940) 

3 Chandrasekhar, S., Proc. Roy. Soc. (London) A, 204, 435 (1951) 

‘ Batchelor, G. K., Proc. Camb. Phil. Soc., 43, 5383 (1947). 

5 Proudman, I. (in press). [ am indebted to Drs. Proudman and Batchelor and to 
Professor Chandrasekhar for having made a manuscript copy of this paper available to 


mec, 


6 Chandrasekhar, S., Proc. Roy. Soc. (London) A, 200, 20 (1949), 


ERRATA 
In the article “On the Origin of the Solar System,”’ these PROCEEDINGS, 
37, 1-14 (1951), the following corrections should be made: 
Page 6, formula (4) for d read a. 


Page &, table 1, last column, for 1.0 read 1.10. 
GERARD P. KUIPER 
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